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Abstract 
The electron-pairing mechanism in unconventional high temperature 
superconductors (HTS) has not been resolved. The author proposed that the 
electron-pairing medium of unconventional HTS is the change of the electron clouds 
of transition metal ions, which is analogous to the lattice vibration in conventional 
superconductors. Real-time evolution of the electron clouds of transition metal ions 
under excitations in La2Fe2zAs2O2, FeSe sheet, Fe2KSe2, CaCuO2, and 
HgBazCazCu30s was calculated by the time-dependent density functional theory 
(TDDFT). The characteristic frequency is about 90-250 meV, which is equivalent to 
the lattice vibration frequencies, showing that the change of the electron clouds of the 


transition metal ions can be the electron-pairing medium in unconventional HTS. 


1. Introduction 


12] and iron-based superconductors [3141 


The discovery of copper-oxide | 
indicates that electron-lattice interactions cannot explain the electron-pairing 
mechanism in unconventional HTS. According to the BCS theory '>!°!, the 
superconducting transition temperature (Te) caused by electron-lattice interactions (at 
normal pressures) cannot be higher than 40K. The pairing mechanism for 
unconventional HTS is still under debate. Mr. Anderson "7! has raised an important 
question in 2007: Is There Glue in Cuprate Superconductors? The author believes that, 


if a medium does exist in unconventional HTS, the medium should change under 


excitations and the characteristic frequency should be close to that of the lattice 


vibrations. Table 1 gives the maximum phonon frequencies and Te of three typical 
conventional superconductors. It is reasonable to guess that the characteristic 
frequency of the electron-pairing medium in unconventional HTS should be about 
100-200 meV. A higher frequency gives a higher Te. If the frequency is too low or too 
high, it is unlikely to be an electron-pairing medium in HTS 


Table 1 Maximum phonon frequencies and T, of three typical conventional superconductors. 


Pp’! MgB"! H,S"! 
œ / meV 9 90 250 


T./K 7 39 164 


Based on the above consideration, the author has studied eight typical 
unconventional superconductors (Fe2KSe2, Laz2Fe2As202, Nd2Fe2As202, Ba2zFe4sAsa, 
YBazCu307, HgBazCa2Cu3Oo, Tl2BazCaCu20s and Bi2Sr2Ca2zCu3O10) l°], Under a 
static electric field, the electron clouds of transition metal ions change significantly. 
The author suggests that the change of the electron clouds of transition metal ions 
may be a new electron-pairing medium. The pairing mechanism is as follows. When a 
free electron comes to a new place, the electron clouds of the neighboring transition 
metal ions will change. In this way, the charge densities around the free electron will 
decrease. When the free electron leaves, the electron clouds of the transition metal 
ions will not relax immediately, so that there will be a region lack of charge. Another 
free electron will be attracted. An attraction appears. This mechanism is essentially 
the same as the electron-lattice interaction, except that the medium is the change of 
the electron clouds, not the displacement of the ions. 

Based on the previous studies, the author has made a further investigation, 
exploring the frequencies of the change of the electron clouds of transition metal ions. 
The real-time evolution of charge densities under some excitation of LazFe2As2O2, 
FeSe sheet, Fez2KSe2, CaCuOz and HgBazCazCu3Os has been calculated by TDDFT 


L12] method. This paper will report the methods and results. 


2. Methods 

Calculations were conducted within the Octopus package [11'S]. The 
generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) was 
used to describe the exchange-correlation energy. The GGA+U method "6ll!7] was 
taking into account to deal with the strong correlation of the 3d electrons. HSCV 
pseudopotentials |'*! were adopted. Approximated Enforced Time-Reversal Symmetry 
(AETRS) algorithm was used to approximate the evolution operator and the time step 


is 0.002 fi/eV. 


3. Results and Discussions 
3.1 La2zFe2As202 

The k-point setting is 4x4x2 and 5p°d!6s* of La, 3p°d°4s* of Fe, 3d!°4s’p? of As 
and 2s’p* of O are taken as valence electrons in the calculation. The change of charge 
densities is induced by a time-dependent electric field: 


fO=EF, cos e io coso If |t- t| >To, then f(t)=0. 


0 
where Fo = 1.0 eV/angstrom, œ = 1.8 eV, to = 5.0/eV and to = 5.0/eV. The excitation 
vanishes after 5000 steps and the charge densities vibrate freely thereafter. After 8000 
steps, the charge densities are recorded. Fig. 1 shows the crystal structure and the 
real-time evolution of the charge density. 100 represents the charge density difference 
of the 8100th step and that of the 8000th step and the same below. The change of the 
electron clouds of the Fe ions become obvious gradually with the evolution steps. 
After 100 steps, almost no change can be seen, and after 500 steps, the change has 
been obvious. The charge density change reaches its maximum value after 800 steps, 
and the corresponding time is 1.6 f/eV. It is the time from zero to the maximum. The 
time of one period should be 4x1.6 A/eV and the frequency is about 160 meV. The 
frequency given in this way is not very accurate and may have an error as large as 
+20%, but the frequency is consistent with the author’s guess, indicating that it can be 


the electron-pairing medium in HTS. 
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Fig.1 Crystal structure and the charge density evolution of La2Fe2As202. The isosurface is 0.02 
num electron/bohr’. The yellow color represents an increase of the charge density, while the blue 
a decrease (the same below). 

3.2 FeSe Sheet [°] 

3p°d°4s* of Fe and 4s’p* of Se are taken as valence electrons. A 4x4x1 k-point 
setting was used and the excitation is the same as above. Fig. 2 shows the results. By 
careful examination, it can be found that the charge density change first increased and 
then decreased, and then increased again. To give a better illustration, Fig. 3 shows 
the maximum value of the charge density change as a function of the evolution steps. 
There are two maximum values: after 550 steps and after 1400 steps. There may be 
two modes in the evolution. This can also be seen from the shapes the charge density 
changes after 600 steps and after 1400 steps. The frequencies are approximately 230 


meV and 90 meV, respectively. 
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Fig.2 Crystal structure and the charge density evolution of FeSe sheet. The isosurface is 0.01. 
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Fig.3 FeSe sheet: The relationship between the maximum value of charge density change and the 
evolution steps 

3.3 KoFes4Se4-2e 

Two valence electrons were subtracted from K2Fe4Se4 and the valence state of K 
is +1, Se is -2, and Fe is +2. The excitation is the same as above. There are also two 
modes in K2Fe4Se4-2e. Fig. 4 shows the results obtained with a k-point setting of 
2x2x1. The maximum of the charge density change is after 500 steps, and the 
frequency is approximately 250 meV. Fig. 5 shows the results with a k-point setting of 
4x4x2. The change reaches its maximum after 1100 steps, and the frequency is 
approximately 100 meV. The difference between the two modes can be seen from the 


patterns of the charge density changes. Different k-point settings may lead to different 
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modes, which needs further investigation. 
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Fig. 4 Crystal structure and evolution of the charge density of K2Fe4Se4-2e. The k-point is set to 


2x2~x1 and the isosurface is 0.04. 
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Fig.5 Crystal structure and evolution of the charge density of K2Fe4Se4-2e. The k-point is set to 


4x4x2 and the isosurface is 0.01. 


3.4 CaCuQ? 7°! 

5pĉd!?6s? of Hg, 5p°6s” of Ba, 3p°4s” of Ca, 3p°d!°4s! of Cu and 2s’p* of O are 
taken as valence electrons. To excite the charge vibration, the position of a Ca atom is 
slightly changed. The original coordinate is (0, 0.5, 0.5), and modified to (0.02, 0.52, 
0.5), ie., the Ca atom was moved in the (001) direction by 0.15 A. Ground state was 
calculated. Then the coordinate was modified back to (0,0.5,0.5) for TDDFT 
calculation. Charge densities were recorded after 13000 steps. Fig. 4 shows the results. 
The electron cloud of the center Cu ion marked by an arrow changes gradually. The 


change reaches its maximum value after1300 steps and the frequency is about 100 


meV. 
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Fig.6 Crystal structure and evolution of the charge density of CaCuOz. The k-point is set to 4x4x2 
and the isosurface is 0.03 

3.5 (HgBa2Ca2Cu308)2 

The k-point is set to 2x21 and vibrations of the charge densities were excited in 
two ways. Fig. 7 shows the results excited by a time dependent electric field. Charge 
densities were recorded after 8000 steps. Except Ca ions, the electron clouds of other 
ions changed obviously. It is worth noting that the oxygen ions in the copper-oxide 
plane undergo obvious changes. The change in the charge densities of almost all the 


ios is maximized after 500 steps with a frequency of approximately 250 meV. 
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Fig. 7 Crystal structure and evolution of the charge density of (HgBazCa2Cu3Os) excited by a time 


dependent electric field. The isosurface is 0.005. 


Fig. 8 shows the results excited by a slight change in the position of an O atom. The 
original coordinate of the O ion in the copper oxide plane is (0.75, 0.25, 0.0), and 
modified to (0.8, 0.3, 0.0), which is equivalent to a change of about 0.4 A along the 
(110) direction. The ground state was calculated. Then modify the coordinate back to 
(0.75, 0.25, 0.0) and conduct TDDFT calculations. After 14000 steps, charge densities 
were recorded. The change of the electron cloud of the Cu ion indicated by the arrow 


reaches its maximum after 500 steps. The frequency is approximately 250 meV. 


Fig.8 Crystal structure and evolution of the charge density of (HgBa2Ca2Cu30s) excited by 


changing an O atom. The isosurface is 0.05. 


4. Conclusion and research prospects 
Real-time evolution of the electron clouds of transition metal ions in 
unconventional HTS was calculated by the TDDFT method. The frequencies of 
evolution of electron clouds are about 90-250 meV, matching fairly well with the 
frequencies of the lattice vibrations at normal and high pressures. Though the 
frequencies obtained by this method are not accurate, it can give a significant 
evidence that the change of electron clouds can be the electron-pairing medium in 
HTS. It is worth noting that there are different modes of the change of the electron 
clouds of the transition metal ions, corresponding to different frequencies. Moreover, 
the electron clouds have very high frequency vibrations, and the author observed the 
vibration with lowest frequency. 
Why do the electron clouds of transition metal ions change this way? The main 


reason is that the 3d shell of transition metal ions is not fully filled. The electron 
9 


clouds possess no spherical symmetry, and easy to change under electric fields. 
Furthermore, for an ion with an incompletely shell, such as O ions in copper-oxide 
planes, same behaviors may occur. 

If the change of the electron clouds of transition metal ions is indeed the pairing 
medium, superconductivity can be divided into two categories according to the 
pairing medium: “displacement of ions” and “electron cloud change of ions”. 

The research of this paper is still incomplete. First, the frequencies should be 
calculated accurately. Second, experimental results must be given; Third, theoretical 


calculation of the Te of unconventional HTS should be achieved. 
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